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Possible golden events for ringdown gravitational waves 
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There is a forbidden region in the parameter space of quasinormal modes of black holes in general 
relativity. Using both inspiral and ringdown phases of gravitational waves from binary black holes, 
we propose two methods to test general relativity. We also evaluate how our methods will work 
when we apply them to Pop III black-hole binaries with typical masses. Adopting the simple mean 
of the estimated range of the event rate, we have the expected rate of 500 yr -1 . Then, the rates 
of events with signal-to-noise ratios greater than 20 and greater than 50 are 32 yr -1 and 2 yr - , 
respectively. Therefore, there is a good chance to confirm (or refute) the Einstein theory in the 
strong gravity region by observing the expected quasinormal modes. 

PACS numbers: 04.30.-w,04.25.-g,04.70.-s 


I. INTRODUCTION 

Black hole (BH) singularities appear unavoidably in general relativity (GR). However, as a physics law the allowance 
of the presence of singularities will not be acceptable even though they are hidden behind the event horizon. Therefore, 
various possibilities of the singularity avoidance have been discussed. Some replacement of singularities is required as 
a complete theory which can describe the BH evolution inside the horizon. Although it is totally unknown how the 
singularities are to be regularized, there are a lot of proposals motivated by the string theory and/or the BH information 
paradox. Some of them, such as gravastars [T], fuzzballs (see, e.g., Ref. [2] for the review) and firewalls 00] change 
the structure of BH spacetime even outside the horizon. Also, an interesting class of singularity and ghost free theories 
of gravity has been proposed by Ref. 00. 

In this paper, we consider binary black hole (BBH) systems, and use gravitational wave (GW) observations as a 
tool to test whether the newly formed black hole genuinely behaves like the one predicted by GR or not. There are 
various methods proposed for testing GR by means of quasinormal mode (QNM) GWs (see an extensive review [7]), 
for example, tests of the no-hair theorem combining two or more modes 0. QNMs dominate the GWs at the ringdown 
phase of BBH mergers (see also Ref. 0). In Ref. ,T(J, testing Hawking’s area theorem IT] has been discussed, which 
is possible if we can determine the masses and spins of BHs before and after merger independently with a sufficiently 
high accuracy. 

One of the methods that we propose in this paper is the following simple one. First, we extract the binary parameters 
of BBHs by taking correlation with the post-Newtonian (PN) templates [T2JDE]- We assume that we know sufficiently 
high PN-order terms to describe the inspiral phase well. Thanks to the development in numerical relativity (NR) [l4l - 
nsi, now we can use simulation results to describe the BBH merger phase, deriving accurate gravitational waveforms. 
Next, if GR is correct, after the merger phase we will observe ringdown (QNM) GWs from the remnant BHs (see, 
e.g., Ref. 0 for a review of the QNMs). If we do not detect the QNMs as expected, it is possible to distinguish the 
remnant object from the BHs that are predicted by GR within the assumptions mentioned above. It should be noted 
that our approach is similar to Ref. nzi. in which the authors discussed the improvement in parameter estimation by 
combining inspiral and ringdown GWs from compact binaries. By contrast, the focus of our work is on the test of 
GR. 

The other method shown in this paper is even simpler. When we focus on the dominant QNM, there is a forbidden 
parameter region in GR. Just using the ringdown GWs, we can directly discuss whether the QNM from the remnant 
compact object is consistent with the one from a BH predicted by GR or not. 

This paper is organized as follows. In Sec. [H] we summarize our tools, the inspiral and ringdown waveforms from 
BBHs, the fitting formulas for the remnant mass and spin, and the matched filtering and parameter estimation in 
the GW data analysis. In Sec. m two simple tests of GR are presented. One is to use only the ringdown GWs, 
and the other is the combination of inspiral and ringdown phases. Finally, we summarize and discuss our approach 
in Sec. |IV| In this paper, we use the geometric unit system, where G = c = 1, and the characteristic scale is 
1 M q = 1.477km = 4.926 x 1CT 6 s. 
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II. PREPARATION 
A. Target of gravitational waves 

According to Kinugawa et al. HME typical total and chirp masses for Pop III BBHs are ~ 6 OM 0 and ~ 30M 0 , 
respectively. Here, the chirp mass of a binary is defined by M. = Mrf / 5 with the total mass M = m\ ± m 2 and 
the symmetric mass ratio 77 = m\m 2 /M 2 . This means that 77 ~ 1/4 for almost equal mass BBHs, which we think 
typical ones. In the following discussion, we focus on equal mass BBHs. Although spins of BBHs can be important, 
we ignore them here for the following reason. If we take into account the spins, one may think that the accuracy of 
parameter estimation might be significantly reduced due to the degeneracy among the orbital parameters. However, 
in that case the orbital precession induced by the spin effects modulates the gravitational waveform. Therefore, to 
a certain extent, this additional information can compensate the loss of accuracy due to the degeneracy. Hence, for 
simplicity, we use only the nonspinning inspiral waveform. 

The inspiral phase of GWs from BBHs has been extensively studied using the PN approximation m ■ If we adopt 
the stationary phase approximation (SPA) [20], we can easily transform the waveform into the expression in the 
frequency domain as A( m e 1 ^ ern . Here, we discuss only the (t = 2, m = 2) mode, and the phase is written as 

«"> = 2 S 1 ' 3 - 2 ®^ I + l2^[ 1 + 0( '' 2 )]' (1) 

where v = (Mirf) 1 / 3 , t c and 4> c are the time and the phase of coalescence, and the higher-order PN terms are 
summarized, e.g., in Eq. (A.21) of Ref. [21] - The appropriate SPA amplitude in the frequency domain is deduced 
from the time domain description A 22 by 

~ 1 7 tM 

A22 = A 22 ^Wv' (2) 

where v is given in Eq. (A. 15) of Ref. [21] . 

After passing the innermost stable circular orbit (ISCO), the BBHs swiftly plunge to merge. Therefore, we terminate 
the inspiral GW analysis at the GW frequency for the (m = 2) mode at ISCO, ,/isco = ( 6 3 // 2 7 r ¥) _1 [221. For a typical 
case with M = 60M 0 , 77 = 1/4, this ISCO frequency is given by /isco = 73.28Hz. 

We can discuss the waveform from the merger phase accurately using NR simulations [T4HT6] . The whole of GW 
waveforms from BBH coalescence are also well modeled in the effective-one-body approach (see, e.g., Ref. [23] for the 
latest development). However, here, we do not make use of the GWs from the merger phase. There is much progress 
in the understanding of the mass, spin and recoil velocity of the remnants after BBH mergers which allows us to 
connect the observation of the inspiral phase to the ringdown phase (see, e.g., Ref. [24] for the latest formulas). Here, 
we use the formulas for initially nonspinning cases. The phenomenological fitting formulas for the remnant mass and 
spin are given by [24] 


M T , 


M 


= ( 477) 2 (M 0 + K 2d 6m 2 + K if 5m 4 ) + 1 + r](Ei SCO + 11) 


Sm 6 , 


Sr, 


^rem — 


M 2 

rem 


= (477) 2 (L 0 + L 2d Sm 2 + L 4f 5m 4 ) + rjJiscoSm 6 , 


(3) 

(4) 


where Sm = (mi — m 2 )/M (= —\J1 — ir] for m\ < m 2 ) and Eisco and Asco are the specific energy and angular 
momentum at ISCO in the test particle approximation (see, e.g., Ref. [23]). Mo, K 2d , K 4 f, L 0 , L 2d and L 4 f are 
the fitting parameters summarized in Table VI of Ref. [24]. a is a/M of the Kerr BH with the mass M and Kerr 
parameter a. More specifically, for equal mass cases, i.e., 77 = 1/4 and 5m = 0, we have 


M T , 


M 


0.951507 ± 0.000030, 
0.686710 ± 0.000039, 


(5) 


including the magnitude of numerical errors. As we noted before, the remnant mass becomes M rem = 57.O9O4M 0 for 
a representative case with M = 60M 0 , 77 = 1/4. 

The above formulas obtained by fitting the results of BBH simulations in the case of nonprecessing BBHs have 
1% relative error, which is mainly caused by the extraction of the GW radiation at a finite radius and finite mesh 
resolution in the NR simulations. The radial extrapolation errors will be reduced by using a perturbative extraction 
method [Ml Efl. Also for precessing BBHs, we may have much larger errors. Although these errors are directly related 
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to the following analysis, we expect that the fitting formulas will be improved by more NR simulations. Therefore, 
we just ignore them in the following analysis. 

Using the estimated remnant BH’s mass and spin, we discuss the ringdown phase. The waveform is modeled as 

fc (t — to) 

— cos (2 Tt fAt- to) - (j>o) for t>t 0 , /gN 

for t <t 0 , 


h(fc, Q, ^0, ^O, t) 


where to and <fo are the initial ringdown time and phase, respectively. The central frequency f c and the quality factor 
Q are related to the real (/r) and imaginary (//) parts of the QNM frequency as 


fR = fc, fl 


A 

2 Q ’ 


(7) 


which depend on the harmonics index (£, m ) and the overtone index n. Here, we focus on the dominant (£ = m = 2) 
least-damped (n = 0) mode and the fitting formulas for f c and Q are given in Ref. [28] as 

fc = 0 * [1-5251 - 1.1568(1 - a rem )° 1292 ] 

Z7T JVl rern ^ 

= 538 ' 4 (^y [1-5251 — 1.1568(1 — a rem )°' 1292 ] [Hz], ( 8 ) 

Q = 0.7000 + 1.4187(1 -a rem ) -0 - 4990 . (9) 


For the fiducial values, M = 60M 0 , 77 = 1/4, we have M rem = 57.O9O4M 0 and a rem = 0.686710, and the above 
formulas derived based on GR predict f c = 299.5Hz and Q = 3.232 for the ringdown GW. Here, it is noted that 
the fitting formulas in Eqs. ^ and Q have 2% and 1% errors, respectively. Therefore, although we use the fitting 
formulas for simplicity in this paper, we should use the original data in Ref. [29] for the strict analysis. 


B. Matched filtering and parameter estimation 


To analyze the GWs from the inspiral and ringdown phases, we use the matched filtering method because the 
waveforms are known well. Using the inner product, 


(a 1 5) = 43? 



5(/)fc*(/) 

Sn(f) 


( 10 ) 


where S n (f) denotes the power spectral density of GW detector’s noise, the optimal signal-to-noise ratio (SNR) for 
a waveform h is given by 


SNR = (/i|h ) 1/2 


= 2 


Hf)f 

SAD 


df 


1/2 


(ii) 


We assume a single GW detector, KAGRA isa [35, here. In Fig. [l] we show the expected noise curve of KAGRA 
[bKAGRA, VRSE(D) configuration] presented in Ref. [32], which can be fit well by 


SAf) 1/2 = 10” 26 (6-5 x 10 10 / -8 + 6 x 10 6 / -2 - 3 + 1-5/ 1 ) [Hz" 1 / 2 ], (12) 

where the frequency / is in units of Hz. Of course, we can discuss the other detectors (Advanced LIGO [35], Advanced 
Virgo [51], GEO-HF [35], and so on) just by changing S n (f)- 

To calculate the parameter estimation errors for the inspiral and ringdown GWs, we use the Fisher information 
matrix, 


r — 


W 


8h\ 

w / 


e=e t , 


(13) 


where 6 l is the parameters of the waveforms and $true denotes the true values of the parameters of the source. Then, 
the rms errors in the estimated parameters and the covariance between two parameters are derived by the inverse 
matrix (T -1 )®! as 

(A Arms = 
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FIG. 1: Fitting curve based on the sensitivity curve of KAGRA [bKAGRA, VRSE(D) configuration] shown in Ref. |32| . 




(14) 


Here, we do not sum over i and j. (A0*)rms scales as 1/SNR. 

For the inspiral phase, we calculate the parameter estimation errors for { M , r], t c , $ c }, Here, we use the total 
mass instead of the chirp mass for the parametrization of the inspiral signal, simply because the fitting formulas 
for the remnant mass and spin are written in terms of M and ry. To evaluate the inner product (101, we take the 
integration range between 10Hz and /iscOj For the ringdown phase, we discuss the parameter estimation with respect 
t° {fc, Q, to, <M, and the frequency interval for the integration is between 10 and 2500Hz. 

We should note that in practice the location of the GW source in the sky and the GW polarization angle in a 
detector frame are also the parameters to describe the GW signals. For example, Ajith and Bose (35] estimated the 
parameter errors of BBHs in a single detector or a detector network for the case of the complete set of parameters. 
This direction to discuss more precise parameter estimation is one of our future studies. 


III. SIMPLE TEST OF GR 


According to Ref. m , individual SNRs for the inspiral and ringdown phase signals are comparable for a gravitational 
wave detector, KAGRA. when the total BBH mass (~ remnant BH mass) is ~ 60 Mg. Since there is a difficulty in 
determining the initial ringdown amplitude due to the ambiguity of the initial time, for simplicity, we set the SNRs 
for the inspiral and ringdown phases to be equal for the typical case (with M = 60 Mq and r] = 1/4 for inspiral 
and M rem = 57.0904Mq and a rem = 0.686710 for ringdown). The assumption of the same SNR for the inspiral and 
ringdown phases is just for simplicity, and we can apply the following analysis for general SNR cases. The information 
of SNRs is imprinted in the Fisher information matrix of each phase. We briefly discuss the effect by setting different 


SNRs for the inspiral and ringdown phases in Sec. IV 


A. Only ringdown 

First, using only the ringdown GWs, we propose a simple method to test whether the compact object emitting the 
ringdown GWs is a BH predicted by GR or not. Figure [2] shows the QNM frequencies for the dominant (t = 2, m = 2) 
least-damped (n = 0) mode in the (/r, fj) plane. In GR, the top-left side of the thick black line is prohibited. The 
boundary thick black line corresponds to the Schwarzschild limit, which is obtained by setting a rem = 0, i.e., 

0.236, (15) 

Jr 

in Eqs. 0. § and ([9]). In principle, if we obtain the parameters in the forbidden region from GW observations, we 
can conclude that the compact object is not the one predicted by GR. 
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FIG. 2: Real (fn) and imaginary (//) parts of QNM frequencies for the dominant (£ = 2, m = 2) least-damped (n = 0) mode. 
The (black) thick line shows the Schwarzschild limit, and the (red) curves are for various mass cases terminated at the spin 
a = 0.998 [38]. From the top of the (red) curves, we are considering BH masses, AI/Mq = 30, 60, 90, 120, 150, and 180, 
respectively. The (red) circles for each line denote the spin dependence a = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 from 
the left. 


However, in practice, there are parameter estimation errors in the GW data analysis. For our typical example with 
f c = 299.5Hz, Q = 3.232, t 0 = 0, and 4>o = 0 [Jr = 299.5Hz and // = —46.34Hz from Eqs. ([7])], we show the contours 
of the parameter estimation errors in Fig. [3] Here, since we do not discuss the errors of t 0 and </>o, we integrated 
the probability distribution over both to and </>o [23J ■ In our typical case, the expected errors are sufficiently small 
to fit the ringdown GW with SNR = 50 within the QNM parameter region allowed in GR at the 5cr level. On the 
other hand, the error circle for the signal with SNR = 20 is not sufficiently small in this sense at that level, while it is 
small enough for 3er level arguments. Here, 5cr (3tr) denotes that for the bidimensional (Rayleigh) distribution, which 
means that the probability falling in the 5er (3cr) circle is about 1 — 3.7 x 10 -6 (1 — 1.1 x 10 -2 ) since the distribution 
has 2 degrees of freedom. [In the case of the ordinary one-dimensional Gaussian distribution, the probability falling 
in the 5<r (3<r) region is about 1 — 5.7 x 10~' (1 — 2.7 x 10 -3 ).] 




f_R [Hz] f_R [Hz] 


FIG. 3: In the (/r, //) plane, the left and right panels show the parameter estimation in the cases with SNR = 20 and 50 for 
the typical case (with A / rem = 57.0904Mq and a re m = 0.686710), respectively. The (black) thick line shows the Schwarzschild 
limit which is same as that in Fig. [2] and the ellipses are the contours of la, 2a, 3 a, 4ct, and 5<r. Here, the time and phase 
parameters (to, <(>o) have been marginalized out. 

To discuss the region prohibited by GR, we present the parameter estimation for the Schwarzschild (a rem = 0) case 
in Fig. [4j Here, we fixed SNR = 50 and considered the remnant masses, M rem /MQ = 45, 60, and 90. From the 5 a 
contours, there is an upper bound of the GR prediction for \fi\/ /r, and we find that the region of \fi\/ /r > F m ax for 
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each mass case is rejected by GR. Here, F max , which denotes the maximum of \fi\/fn allowed in GR, is 0.321 (for 
Ikf rem = 45 Mq), 0.320 (60M 0 ) and 0.316 (90M Q ) for SNR = 50. If NR simulations for the extreme spinning BBH 
are available, we can also give the lower bound of the GR prediction for |//|//r. 



FIG. 4: In the (/r, //) plane, this figure shows the parameter estimation in the cases with SNR = 50 for a Schwarzschild black 
hole with M rem /Mg = 45 (right), 60 (center) and 90 (left). The (black) thick line shows the Schwarzschild limit which is same 
as that in Fig. [2] and the ellipses are the contours of la, 2a, 3ct, 4ct, and 5cr. Here, the time and phase parameters (to, 4>o) have 
been marginalized out. 

It is noted that a powerful method to find ringdown signals in multiple GW detectors has been proposed by Talukder, 
Bose, Caudill and Baker flU]. Although we have considered the above GW data analysis with a single detector, we 
may expect a better parameter estimation in a detector network. 


B. Consistency analysis with inspiral and ringdown 


Next, we propose a consistency test by combining the data from inspiral and ringdown GWs. We use the PN 
waveform for the inspiral phase to extract the binary parameters, and the formulas in Eqs. ([3]) and Q of Sec. [TT] are 
applied to obtain the GR prediction for the parameters of the remnant black hole. Then, we can present the QNM 
frequency expected in GR in the (/r, //) parameter space. 

To take into account the observational errors in the estimate of the expected QNM, we assume that the true signal 
is given by the GR template with 6, and the parameters estimated from the inspiral and ringdown signals are 0i nS p 
and 0Rmg, respectively. Here, 6 consists of the parameters {f c , Q, to, <j>o}, which are commonly used for the ringdown 
GW data analysis. For the ringdown phase, we treat the above parameters to calculate the parameter estimation 
errors and assume the Gaussian distribution for the parameters. In the inspiral-phase analysis, we use another set of 
parameters 6 = {M, 77 , t c , $ c }. 

Here, it is useful to have the relation between the inspiral parameters 0 and the ringdown parameters 6 as fitting 
functions. From Eq. Q, we have 

“~ m = a880028 (5^5) - 0 148761 + 000180881 (s(hf ■ (16) 


The above relation gives one-to-one mapping in the parameter ranges, 0 < 77 < 0.35282872 and 0 < a rem < 0.99800367. 
It is noted that, although 77 > 0.25 is an unphysical value, we allow the values here. Combining the above equation 
with Eq. ([9]), we find that 77 is fitted as a function of Q to obtain 


1 


0.353039 


0.208266 


10.9583 21.4540 

—Q* W~ 


(17) 


The restriction on the parameter space to keep the one-to-one mapping becomes 2.11870 < Q < 32.2555. The decay 
time is calculated as Q/(nf c ). Using Eqs. <§ and Q (and also the above fitting functions for a lem and 77 ), the total 
mass M in the inspiral phase is written by f c and Q as 


1 


f 1 \ 

0.06461671 

Tc 

—0.0434932 - 0.127430 In | 

v q +0 - 163772 J 

VQ 


( 18 ) 
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To find the expected parameter region of the QNM, we use the following simple estimator (more detailed studies, 
e.g., by using Markov chain Monte Carlo methods, will be presented in future): 


F(9) =Af exp 


_i f insp ( ^( 0 ) _ r(e lnsp )W(e) - oHe i nsp )) - *r* ing (^ - <4 inK )(<p - ^ Rine ) 


2 *■* 


Ring/ 


(19) 


where AT is a normalization constant which we do not take care of, and F*” sp and r^ mg denote the respective Fisher 
information matrices after integrating the probability distribution over (t c , < f> c ) and (to, 

The strategy to estimate the parameter region by using Eq. (19) is as follows: 


(1) For given 0(0 Insp )(=: $i nsp ) (in practice, we give 0i nsp = {M = 60M 0 , r) = 1/4} and derive 0i nsp ), we calculate 


fj" sp with the bKAGRA noise curve. 

L J 


(2) Assuming the narrow ringdown signal in the frequency domain, we prepare r^ mg for the white noise (analyti¬ 
cally). 


(3) For given 0ru ng (and F^ ms for it), we find the maximum of Eq. (19) by 


dF{9) 

dd l 


= 0 . 


( 20 ) 


(4) Inserting the solution of the above equation 0 = {/ c , <5} back into Eq. (19), we check whether the situation 


with the parameters (0i nsp , 9 Ri ng , 9) is in the 5cr level of the detector noise realization or not. 


Here, 5er denotes that the value in the exponent in Eq. (19) becomes —5 2 /2, which means that the probability 
falling in the 5<r circle is about 1 — 5 x 10~ 5 since the distribution has 4 degrees of freedom. Employing our fiducial 
values, M = 60Afg, q = 1/4, the expected region of the QNM frequency in the 5cr level is shown in Fig. [5 Here, we 
have fixed SNR = 50 for both the inspiral and ringdown GWs. Compared with the right figure of Fig. [3j the allowed 
region in this figure has a larger extension in the horizontal direction. This is due to the parameter estimation errors 
of the inspiral phase. We repeat the meaning of this plot. Under the condition that we measure the values of both 
9i nsp and 0Ri ng , we choose the most probable values for the parameters 9. Assuming that the true values are the 
most probable values as used in the usual Fisher-matrix analysis, we evaluate the probability that the detector noise 
produces such a deviation in the measurement of f?i nsp and 9 Ri ng . The probability that the noise realization falls 
outside the contour is 5 x 10 -5 . Therefore, if we find that the parameter estimate from the ringdown signal deviates 
from the prediction from the inspiral signal exceeding the contour in Fig. [3j we can conclude that there is something 
wrong with the GR prediction. Here, under an assumption that the nonlinearity of GR is correct for the inspiral and 
merger phases, it is possible to distinguish the remnant object from the BHs that are predicted by GR. 



FIG. 5: The QNM frequency expected in GR (5<r level) from the inspiral phase with the total mass M = 60Mq and symmetric 
mass ratio r/ = 1/4 [the (red) filled region]. The (black) thick line shows the Schwarzschild limit which is same as that in Fig. [2] 











IV. SUMMARY AND DISCUSSION 


In this paper, we mainly focused on a specific BBH with the total mass M = 6 OAf 0 and the symmetric mass ratio 
77 = 1/4, which would be the typical one for Pop III BBHs T 8 S [19]. It is found that we can perform meaningful tests 
of GR, assuming that the GW signal has SNR = 50. An easy extension of the present study is to treat various total 
mass cases. For total masses lower than M = 60 Mq, we have fewer SNRs for the ringdown phase than those for the 
inspiral phase and expect that a larger elongation in the vertical direction in the (Jr, //) plane because of Fig. [3j On 
the other hand, for total masses higher than M = 60M 0 , we will have a larger elongation in the horizontal direction. 
We also need to discuss various mass ratios and spins in the inspiral phase. The statistical treatment will be also 
improved in our future work. 

In Fig. [5] we have observed that the expected region shows a large elongation in the horizontal direction. This is 
due to the parameter estimation errors for the inspiral signal and, more specifically, originates from marginalizing t c 
and "Fc in the probability distribution. The parameter estimation errors of t c and Tc arise from the short frequency 
integration interval between 10Hz and /isco = 73.28Hz . The number of GW cycles during this frequency range 
is -/Vqw ~ 30. When we change the lower integration bound to 20 Hz, the situation becomes much worse, i.e., the 
number of GW cycles is just IVqw « 6 . 

Here, if we can also detect the inspiral phase by using a space-based GW detector, such as DECIGO BII. the 
situation will improve a lot (see, e.g., Ref. |42| for the synergy in the parameter estimation of binary inspirals). For 
example, Nq w ~ 5400 from 0.5Hz in our specific case. Therefore, even if we assume the same SNR for the inspiral 
phase, the parameter estimation of M and 77 and the QNM prediction will be very precise. 

Kinugawa et al. [15] showed that the expected detection rate of BH-BH mergers by KAGRA with typical total mass 
~ 60 Mq is given by 


262eventsyr 1 2 (SFRp/(10 2 5 M 0 yr 1 Mpc 3 4 5 6 7 )) • Err sys , (21) 

where SFR p and Err sys are the peak value of the Pop III star formation rate and the systematic error with Err sys = 1 
for their fiducial model, respectively. They have estimated that Err sys ranges from 0.056 to 2.3 due to the unknown 
parameters such as the common envelope parameter, the kick velocity, and the loss fraction as well as the unknown 
distribution functions such as the initial mass function and the initial eccentricity function. The minimum value 
corresponds to the worst model in which they adopt the most pessimistic values of the parameters and distribution 
functions within the ranges that are likely. The factor (SFRp/(10 -25 M 0 yr -1 Mpc -3 )) also depends on the models 
and Kinugawa et al. [19] argued that it ranges from 0.019 to 16. Therefore, the event rate of Pop III BH-BH mergers 
which will be detected by KAGRA ranges from 0.28 to 9641 yr -1 . The event rate for Advanced LIGO and Advanced 
Virgo will be similar. Since no such event has been found so far, the event rate should be smaller than 1000 yr -1 . 
Adopting a simple geometric mean of this allowed range, we have a rough estimate of the expected rate of 500 yr -1 . 
Then, the rates of events with SNR > 20 and SNR > 50 are 32 and 2 yr -1 , respectively. Therefore, there is a good 
chance to confirm (or refute) the Einstein theory in the strong gravity regime by observing the expected QNMs. 


Acknowledgments 

This work was supported by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) Grant- 
in-Aid for Scientific Research on Innovative Areas, “New Developments in Astrophysics Through Multi-Messenger 
Observations of Gravitational Wave Sources”, No. 24103006 (H.N., T.T., T.N.) and by the Grant-in-Aid from MEXT 
of Japan No. 15H02087 (T.T., T.N.). We gratefully acknowledge all participants in “Gravitational Wave Physics and 
Astronomy Workshop (GWPAW) 2015,” held June 17 20, 2015 in Osaka, Japan. H. N. would like to thank Y. Nishino 
for useful suggestions. 


[1] P. O. Mazur and E. Mottola, gr-qc/0109035 

[2] S. D. Mathur, Fortsch. Phys. 53, 793 (2005) |hep-th/0502050 . 

[3] S. L. Braunstein, S. Pirandola, and K, Zyczkowski, Phys. Rev. Lett. 110, 101301 (2013) arXiv:0907.1190 [quant-ph]]. 

[4] A. Almheiri, D. Marolf, J. Polchinski, and J. Sully, JHEP 1302, 062 (2013) |arXiv:1207.3123 [hep-th]]. 

[5] L. Modesto, Phys. Rev. D 86, 044005 (2012) arXiv:1107.2403 [hep-th]]. 

[6] T. Biswas, E. Gerwick, T. Koivisto, and A. Mazumdar, Phys. Rev. Lett. 108, 031101 (2012) arXiv:1110.5249 [gr-qc]]. 

[7] E. Berti et al ., arXiv:1501.07274 [gr-qc]. 



9 


[8] O. Dreyer, B. J. Kelly, B. Krishnan, L. S. Finn, D. Garrison, and R. Lopez-Aleman, Class. Quant. Grav. 21, 787 (2004) 
gr-qc/0309007 . 

[9] E. Berti, V. Cardoso, and A. O. Starinets, Class. Quant. Grav. 26, 163001 (2009) arXiv:0905.2975 [gr-qc]]. 

[10] S. A. Hughes and K. Menou, Astrophys. J. 623, 689 (2005) astro-ph/0410148 . 

[11] S. W. Hawking, Phys. Rev. Lett. 26, 1344 (1971). 

[12] L. Blanchet, Living Rev. Rel. 17, 2 (2014) arXiv: 1310.1528 [gr-qc]]. 

[13] G. Schaefer, Fundam. Theor. Phys. 162, 167 (2011) arXiv:0910.2857 [gr-qc]]. 

[14] F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005) gr-qc/0507014 . 

[15] M. Campanelli, C. O. Lousto, P. Marronetti, and Y. Zlochower, Phys. Rev. Lett. 96, 111101 (2006) gr-qc/0511048 . 

[16] J. G. Baker, J. Centrella, D. I. Choi, M. Koppitz, and J. van Meter, Phys. Rev. Lett. 96, 111102 (2006) gr-qc/0511103 . 

[17] M. Luna and A. M. Sintes, Class. Quant. Grav. 23, 3763 (2006) gr-qc/0601072 . 

[18] T. Kinugawa, K. Inayoshi, K. Hotokezaka, D. Nakauchi, and T. Nakamura, Mon. Not. Roy. Astron. Soc. 442, 2963 (2014) 
arXiv: 1402.6672 [astro-ph.HE]]. 

[19] T. Kinugawa, A. Miyamoto, N. Kanda, and T. Nakamura, arXiv:1505.06962 [astro-ph.SR]]. 

[20] T. Damour, B. R. Iyer, and B. S. Sathyaprakash, Phys. Rev. D 63, 044023 (2001) [Erratum-ibid. D 72, 029902 (2005)] 
gr-qc/0010009 . 

[21] P. Ajith, M. Boyle, D. A. Brown, S. Fairhurst, M. Hannam, I. Hinder, S. Husa, B. Krishnan, R. A. Mercer, F. Ohme, 
C. D. Ott, J. S. Read, L. Santamaria, and J. T. Whelan, arXiv:0709.0093 [gr-qc]. 

[22] C. Cutler and E. E. Flanagan, Phys. Rev. D 49, 2658 (1994) gr-qc/9402014 . 

[23] A. Taracchini et al., Phys. Rev. D 89, 061502 (2014) arXiv:1311.2544 [gr-qc]]. 

[24] J. Healy, C. O. Lousto, and Y. Zlochower, Phys. Rev. D 90, 104004 (2014) arXiv: 1406.7295 [gr-qc]]. 

[25] A. Ori and K. S. Thorne, Phys. Rev. D 62, 124022 (2000) gr-qc/0003032 . 

[26] H. Nakano, Class. Quant. Grav. 32, 177002 (2015) arXiv: 1501.02890 [gr-qc]]. 

[27] H. Nakano, J. Healy, C. O. Lousto, and Y. Zlochower, Phys. Rev. D 91, 104022 (2015) arXiv:1503.00718 [gr-qc]]. 

[28] E. Berti, V. Cardoso, and C. M. Will, Phys. Rev. D 73, 064030 (2006) gr-qc/0512160 . 

[29] http: //www.phy. olemiss . edu/~berti/qnms .html 

[30] K. Somiya (KAGRA Collaboration), Class. Quant. Grav. 29, 124007 (2012) arXiv:1111.7185 [gr-qc]]. 

[31] Y. Aso et al. (KAGRA Collaboration), Phys. Rev. D 88, 043007 (2013) arXiv: 1306.6747 [gr-qc]]. 

[32] http://gwcenter.icrr.u-tokyo.ac.jp/researcher/parameters 

[33] J. Aasi et al. (LIGO Scientific Collaboration), Class. Quant. Grav. 32, 074001 (2015) arXiv: 1411.4547 [gr-qc]]. 

[34] F. Acernese et al. (VIRGO Collaboration), Class. Quant. Grav. 32, 024001 (2015) arXiv: 1408.3978 [gr-qc]]. 

[35] C. Affeldt et al., Class. Quant. Grav. 31, 224002 (2014). 

[36] P. Ajith and S. Bose, Phys. Rev. D 79, 084032 (2009) |arXiv:0901.4936 [gr-qc]]. 

[37] N. Kanda (LCGT Collaboration), arXiv:1112.3092 [astro-ph.IM]. 

[38] K. S. Thorne, Astrophys. J. 191, 507 (1974). 

[39] L. S. Finn, Phys. Rev. D 46, 5236 (1992) gr-qc/9209010 . 

[40] D. Talukder, S. Bose, S. Caudill, and P. T. Baker, Phys. Rev. D 88, 122002 (2013) arXiv:1310.2341 [gr-qc]]. 

[41] N. Seto, S. Kawamura, and T. Nakamura, Phys. Rev. Lett. 87, 221103 (2001) |astro-ph/0108011 . 

[42] R. Nair, S. Jhingan, and T. Tanaka, arXiv: 1504.04108 [gr-qc]. 






